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Abstract 

By applying a unitary transformation method, we have derived the leading-order corrections on 
the effective Hamiltonian of a dynamical model developed in Phys. Rev. C54, 2660 (1996) for 
electromagnetic pion production reactions. The resulting energy-independent one-loop corrections 
on the baryon masses and the "fN — > A vertex interaction are associated with the structure of the 
nucleon and A and have been calculated within a constituent quark model. We find that the 
one-loop corrections on the magnetic Ml transition of the "fN — > A are very small, while their 
contributions to the electric E2 and Coulomb C2 transitions are found to be in opposite signs of 
that due to pion cloud effects associated with the scattering states. Our results further indicate that 
the determination of the nonspherical L = 2 components of the constituent quark wavefunctions 
of N and A from the extracted empirical E2 and C2 form factors requires a rigorous and complete 
calculation of meson cloud effects. We also find that the one- loop corrections on the non-resonant 
pion production operator can resolve the difficulty in describing the near threshold jp — * tt°p 
reaction. Possible future developments are discussed. 

PACS numbers: 13.40. Gp, 13.60.Le, 14.20. Gk, 24.10-i 
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I. INTRODUCTION 



In the past few years, extensive and precise data of electromagnetic meson production 
reactions have become available and some of these data have been used to extract the 
information about the nucleon resonances On the other hand, theoretical models for 
analyzing these reactions are still far from complete. Even in the simplest and well-studied 
A excitation region, none of the most often applied models 0, 0, U, 0, Q has been able to 
give predictions which agree perfectly with the single pion production data accumulated 
recently, in particular the data on spin observables and longitudinal-transverse interference 
cross sections. While these models can give an overall good description of fairly extensive 
data, efforts must be made to remove the remaining discrepancies such that a complete 
understanding of the A resonance can be obtained. The experiences gained from these 
efforts will undoubtly be very useful for investigating the much more complex higher mass 
iV* resonances. In this work, we report on the progress we have made in this direction, 
focusing on the dynamical model we have developed in Refs-P,!^ (called the Sato-Lee (SL) 
model in the literatures). In particular, we would like to explore how the bare jN — > A 
parameters extracted within the SL model can be better understood in terms of the structure 
of N and A. We would also like to see how the non-resonant pion production operator in 
the SL model can be improved. 

We first recall one of the most interesting results from the SL model. It was found that the 
pion cloud effects give very large contributions to the 'yN — > A transition form factors and 
is the source of the differences between the values predicted by the conventional constituent 
quark model and that extracted from empirical amplitude analyses. The predicted very 
pronounced Q 2 — dependence in electric E2 and Coulomb C2 transitions have motivated 
several recent experimental efforts. These pion cloud effects are calculated from the following 
expression 

IW W «) = rVM + / dkk^ {q , k) w _ K(k) 1 _ EN(k) + . r, N Mw, k) 

where rj^y^ is the bare vertex, v^ e is the non-resonant — > ttN amplitude calculated 
from the standard Pseudo- Vector Born terms and the p and ui exchanges, and T vN ^{W, k) 
is the dressed A — > nN vertex. One observes from the above equation that these pion cloud 
effects are due to pions in the scattering states which can reach the on-shell momentum 
asymptotically. 



We now examine how the above procedure is related to our current understanding of 
hadron structure. Because of the chiral symmetry of QCD is spontaneously broken, it is 
generally believed that in the region where the momentum transfer is not too large the 
structure of the nucleon and A can be considered as systems made of constituent quarks 
and virtual pions. We thus expect that their responses to the external electromagnetic 
field can be from the constituent quarks and also from the virtual pions. Obviously the 
pion-loop integration in the above equation do not account for all of the effects due to the 
virtual pions in hadrons. The leading term T® NA must still contain some effects due to 
virtual pions which never go on-shell during the N-A transitions. In this work, we will 
show how the corrections due to these virtual pion cloud effects can be derived by applying 
the unitary transformation method. In a consistent derivation, the one-loop corrections on 
the non-resonant pion production operator of the SL model have also been derived. These 
one-loop corrections are also energy-independent and are different from those due to pions 
in scattering states. These corrections are expected to have important effects in the region 
where the pion electromagnetic reactions are sensitive to the non-resonant amplitudes. 

In section II, we recall a dynamical formulation within which the leading order one-loop 
corrections on the effective Hamiltonian of the SL model are derived. In section III, the 
consequences of these leading order corrections on the 7 A — > A transitions are calculated 
and interpreted within a constituent quark model. The one-loop corrections on the non- 
resonant pion production operator are then investigated in section IV, focusing on the s-wave 
amplitude of the near threshold tt° photoproduction reaction. Possible future developments 
are discussed in section V. 



II. FORMULATION 

As explained in Ref.Q, the SL model is constructed by applying a unitary transforma- 
tion method to deduce from relativistic quantum field theory an effective Hamiltonian for 
describing meson-barvon reactions. The details of the employed unitary transformation has 
been given in Refs. [2|, [7[ and will not be repeated in this paper. Here we only emphasize 
that the starting point of the unitary transformation method is a field theoretical Lagrangian 
density. This is identical to other more familiar approaches for constructing dynamical mod- 
els of meson-baryon interactions, such as those based on the ladder Bethe-Salpeterj^, 0] or 
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three-dimensional ladder Bethe-Salpeter equations |l(3, 11,13, 13|, |14(. In the lowest order, 
all approaches yield very similar, if not completely identical, scattering amplitudes. Their 
differences are in the resulting dynamical equations which are used to include nonperturba- 
tively certain classes of higher order effects that are deemed to be important for the processes 
considered. 

To illustrate the unitary transformation method, it is sufficient to consider a model La- 
grangian density L(ipN,ipA, (fin) describing the pseudo- vector coupling between ir, N and 
A fields. By using the standard canonical quantization procedure, a Hamiltonian can be 
constructed. To simplify the presentation, the spin and isospin variables as well as the anti- 
particle components are suppressed here. The resulting Hamiltonian can then be schemati- 
cally written as 

H = H Q + Hj + H em , (1) 

with 

H = E / dpb B {p)b B {p)E B {p) + f dkal(k) a7T (k)E n (k) 

where b B {p){b B {p)) is the creation(annihilation) operator for a baryon with momentum p, 
and (k)(a(k)) for a pion with momentum k. The energy is defined as E a (p) = (rria+p 2 ) 1 ^ 2 
with m a denoting the mass of particle a. Clearly, H is the sum of free energy operators for 
baryons(I? = N, A) and pion(7r). The strong interaction Hamiltonian in Eq. (pj is 

Hi = £ K B , B + h.c] , (2) 

B,B' 

with 

t Ib',b = J dpdkb B ,(p - k)bB{p)al(k)F nB ',B{p - k, k; p), (3) 

where F^b 1 , B {p'k\ p) is a vertex function describing the strength of the nB <-» B' transition 
illustrated in Fig. [T] The corresponding electromagnetic interaction deduced from applying 
the minimum substitution on the considered pseudo-vector coupling Lagrangian density 
L(i(:n, ^a, 0tt) can be written as H em = J dxA ■ J, where J M is the current density operator 
and A is the photon field. The resulting electromagnetic current can be schematically written 
as 

_ jm _|_ J^, B 4- J^ / B 7T , (4) 
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where J%, J B , B , and J B , B n define •yn — > tt, yB — > B', and the contact jB — > ttB' transitions 
respectively, as illustrated in Fig. |21 The details of these currents will be given in section 
III. 
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FIG. 1: The vertex interaction T^ B , B 
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FIG. 2: Electromagnetic interaction H em : (a) J^ 4 , (b) J B , B , and (c) J B , B . 



The first step of the derivation is to decompose the strong interaction Hamiltonian into 
two terms 



H T = iff + H? 



(5) 



where 



h[ = rl NA + h.c. , (6) 

h? = K N>N + r° A)JV + r° AiA ] + h.c. . (?) 

Obviously, Hf describes the physical process, while the processes in H® can not occur in 
free space because of the violation of energy conservation. The second step is to perform 
unitary transformations on H to construct an effective Hamiltonian, which does not contain 
unphysical processes such as those due to H® . Keeping only the terms up the second order 
in Hi, the resulting effective Hamiltonian is of the following form 



H eff = uIu\hu x u 2 



H + Hf + H? + [UMHemUM + A Q 



(8) 



where U n = exp(iS n ) is the n-th unitary transformation with S n oc (Hi) n , H p has been 
defined in Eq. © and 

H? = ([H?,iS 1 ] + ~[H?,iS 1 ]) p (9) 

Note that the commutators in Eq. © can generate both physical and unphysical processes 
and only the terms for physical processes are kept in H P . The unphysical processes in H e ff 
is contained in the last term of Eq. (JHJ) 

A Q = {[H , iS^ + H?} + {[H , iS 2 ] + H?} , (10) 

with 

H? = {[H?,iS 1 ) + ±[H?,iS 1 ])<>. (ii) 

H® is defined by the same commutators in H P except that only the unphysical processes 
are kept here. 

The desired effective Hamiltonian is obtained by eliminating the unphysical processes 
Eq. (jHJ). Obviously, this can be achieved by imposing the following conditions 

[#o,i£i]+i?? =0, (12) 
[Ho,iS 2 ]+H® =0. (13) 

To find Si, consider the matrix elements of Eq. (j!2)l between any two eigenstates | a > and 
| b > of Hq] for example H \ N >= E N \ N > and H | irN >= {E„ + E N ) \ tcN >. We 
then obtain a relation (E b — E a ) < a \ iSi \ b >=< a \ H® \ b >, indicating that 5*i plays 
the same role as Hf in defining the interaction mechanisms. It is then easy to verify that 
the general solution of Eq. (|T2*|) can be written as the following operator form 

S 1 = -i I dpdk f^'zf^L k ^\ P \ (u J {m^ + m B , - m B )b ] B ,{p - k)b B {p)al{k) 
Z7~LJ Er(v) — En>\V - k) - EJk) 



B'B ' 



+ h.c. (14) 

where the step function is defined as 9(x) = 1(0), for x > (<)0 

For investigating ttN scattering and pion photo- and electro-production at energies below 
two-pion production threshold, it is sufficient to consider interactions defined within the 
Hilbert space iV© A©7riV©7iV. By using Eq. and Eq. (JHJ), we can evaluate the matrix 



elements of , defined by Eq. between two one-baryon states. This will generate 
the one- loop corrections, TP N and £ A , to the masses of N and A, as illustrated in Fig. 
Explicitly, we find 



yO 

N 



yO 
^A 



- J2 [<N\T°J B ^ N \nB'xnB'\S l \N> 

Z B'=N,A 

< N I S 1 I tiB' >< tiB' I T° B I N >], 



2[< A l r °A,Al^X7rA|5 1 |A> 
< A | Si | ttA >< ttA I T° A A I A >]. 



(15) 



(16) 



Using the solution Eq. (JT3J) for Si to evaluate the above two equations, we will get expres- 
sions involving one-loop integrations over energy — independent propagators which are also 
specified in Eq. (|14j) . The detailed forms will be given in the next section where we will 
perform calculations using a model for the vertex interaction T^ B , B . Note that £ A does 
not include the loop over intermediate ttN state since the effects due to A — * ttN is already 
accounted for by if-f of Eq. © and must be excluded in Q interactions. 
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FIG. 3: One-loop corrections S^r and S A on the nucleon and A 
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FIG. 4: 7riV interactions. 

Taking the expectation value of between two ttN states, we then generate the ttN 
potential v n N, illustrated in in Fig. HJ Extending the procedure described above to also 
include spin and isospin indices as well as the anti-particle components and p meson, the 
matrix elements of v n N given explicitly in the SL model can then be obtained. On the other 
hand, the one-loop corrections and S A are not treated explicitly in SL model. 

To determine 5*2 from Eq. (|13|). we need to know the mechanisms contained in defined 
by Eq. (fTTj) . With the solution Eq. (fT4"j) for Si, one can easily see that can generate the 
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FIG. 5: Unphysical processes due to H9. 



unphysical B — > -ktxB processes illustrated in Fig. |5J With the similar procedure employed 
in solving Eq. (fT2*|) for Si, we find that the solution of Eq. (jlHjl from eliminating the 
unphysical processes illustrated in Fig. |S] can be written explicitly as the following operator 
form 

Q • V- I a jujm F "B',b(p -k-k', k'; p - k)F nB B (p - fc, fc; p) 
^2 = -« L / dpdkdk - — — — — ^m7\ f rut\ 8{ 2m n + m B' ~ m B ) 



B , BB ^ E B (p) -E B ,(p-k- fc') - K(fe) - 

9(m n + m B i — m B ) 9(m n + m B — m B 



■E B ,(p-k-k')-E s (p-k) + E„(k'y 2 
9(171^ + m B - m B ) ,0(m n + m B > - m B ) 



+ 9(-m 7T - m B + m s )) 



E B (p)-E B (p-k)-E„(ky 2 
&{,,(p - fe - fe / )6 fl (p)4(fe')ai(fe) + h.c. 



+ 9(-m 7T - m B > + 7715))] 



(17) 



We note that S*2 does not play any role in generating effective the Hamiltonian up to the 
second order in T^ B B , . But it is needed to evaluate the effective electromagnetic interaction 
operator defined by the term \UlUiH ern UiU2\ in Eq. (jHJ. 
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FIG. 6: Leading order terms of pion photoproduction: r° JV)A =(a), v*™ e = (b) + (c) + (d) + (e) 

Keeping only the terms up to the same second order in Hj, we can write [L^^f-f^m^i^] : 
/ dxA ■ J e ff with the effective current defined by 



r eff = j» + [j»,iSi] + -[[j> i ,iSi},iSi\ + [j> i ,is 2 } + [[j»,iSi},is 2 }. 



(18) 



By using the properties of S\ and the electromagnetic coupling illustrated in Fig. one can 
see that the first two terms of Eq. (JT8|l for J^j generate the tree mechanisms shown in Fig. 
EJ Explicitly, we can see the following correspondences : 

Fig. 6a : < A | J^ N | 7 iV >, (19) 

(Fig.6b + Fig.6c) : < nN \ [J^ N , iS t ] + JfaiSi \ 7 iV >, (20) 

Fig.6d : < 7tN I [J^iS x \ | 7 iV >, (21) 

Fig.6e : <nN\ | 7 iV > . (22) 

Extending the procedure described above to also include spin and isospin indices as well 
as the anti-particle components and p and u meson-exchange, the matrix elements for v 1 n 
given explicitly in the SL model can then be obtained. 

The one-loop corrections on the 7 5 — > B' vertex and non-resonant 7 iV — > ttN amplitude 
can be generated from the following operators in Eq. (|18|). 

jKi-ioop = [jm lSi] + 1 [[jm iSl ] t iSl] + [jm z52 ] + [[jm iS y . (23) 

For 7 iV — > A, the possible intermediate states involved in evaluating the one-loop corrections 
are illustrated in Fig. [7|with the following correspondences: 

[J^, B ^iS 1 ]\ 1 N>, (24) 
-iSxJ&^iSx | 7 iV >, (25) 
^[[J^iS^iS^ + iJ^iSi] | 7 iV>, (26) 

^(J^^i5ii5 , i + iSJSxJ^B) | 7 JV > . (27) 

Similar expressions and diagrams are also for the one- loop corrections, < iV | J^~ loop \ 7 iV >, 
for the nucleon electromagnetic form factors. 

The one-loop corrections on the non-resonant 7 iV — > ttN amplitudes can also be obtained 
by taking the matrix element of Jj l ~ loop between nN and 7 iV states. We will elaborate this 
more complex object in section IV. 

With the above derivations, the effective Hamiltonian Eq. (jSJ) within the subspace iV © 
A © nN © 7 iV can be written as 

H eff = [H + Z° N + E° A ] + [T nNA + T jNA ] + [v nN + v in ]. (28) 
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FIG. 7: Loop Correction r^'^ p on 7./V — » A transition. 

The mass correction terms, and S^, are illustrated in Fig. El The vertex interactions 
in Eq. (J2HD are 



tN,A — J- 7r/V,A' 
7 AT, A 



I - p0 1 pi— ioop 

1-vAT.A — J- -yjV,A "T" 7JV,A • 



(29) 
(30) 



Here we have defined 



F° 

1 7 AT, A 



J dx < A \ A - J£ N I > . 



(31) 



The one-loop corrections r^^ p are defined by Eqs. (|23)l - (j27jl and illustrated in Fig. [7| 
Note that up to the second order in Hj, there is no one-loop correction to the ttN — > A 
vertex in Eq. (|29|) . 

The 7riV potential in Eq. (|2*K|) is illustrated in Fig. 0] The non-resonant 7^ — ► 7riV 
transition interaction is defined by 



„7T _ 7 ,*ree 1 l-loop 



(32) 



where f*^ ee is defined by Eqs. (|19| ) -()22| 1 and illustrated in Figs. (6b)- (6d), and t>*~' oop is 
the one-loop corrections which can be calculated by taking the matrix element of Eq. 
between ^yN and ttN states. 
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The SL model can be obtained from the effective Hamiltonian H e ff of Eq. ()28j) by making 
the following simplifications. First, the mass correction terms and £ A are not treated 
explicitly and are included in the physical nucleon mass 771^=938.5 MeV and wia = 1299 
MeV determined in Ref.j^j]. Second, the one-loop corrections Tl~^°£ p are not calculated 
explicitly and r 7 jv,A; instead of T^ ^, is adjusted to fit the data. Finally, the non-resonant 
vl i ~ loop is neglected. 

The above derivation indicates that the method of unitary transformation has provided 
a systematic way to improve the SL model. In the next two sections,we will explore the 
consequences of these leading-order corrections derived in this section. 



III. ONE-LOOP CORRECTIONS ON THE ONE-BARYON PROCESSES 

To evaluate the one- loop corrections Eqs. ()15j) - (fT7)J) for the baryon masses and Eqs. ()24j) - 
()27|) for the 7A — > N and 7A — * A transitions, we need to define the vertex function 
F-kB',b of Eq. (JHJ) and the matrix elements of currents of Eq. (j3J). As an exploratory step, 
we assume that these can be calculated from a model within which the pion is coupled to 
constituent quarks by the usual pseudo-vector coupling and the electromagnetic interaction 
is introduced by the minimum substitution. We further assume that the constituent quarks 
in TV and A are nonrelativistic and only have L = s-wave configurations. Accordingly, the 
usual nonrelativistic limit is also taken to define the couplings of tc and 7 with constituent 
quarks. With these simplifications, we can cast the resulting irB — > B' vertex into the 
following form 

F^ B , BI {p,k ]P ') = 1 1 UbB ' {S B ,, B -k){T B , )B .li)F B ,, B {k) (33) 

Here J* is a vector associated with the pion isospin state i, and F B ',B{k) is a form factor 
calculated from quark wave functions. The spin and isospin operators Sb>b,Tb>b are de- 
fined as follows. For diagonal spin operators they are twice of the spin angular momentum 
operator. 

S NN = 2J = a, (34) 
Saa = 2 J = Sa', (35) 
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while the transition spin operators are defined as 



San = S, (36) 
S NA = S+. (37) 

Within the considered SU(6) quark model, these operators are related to each other , as 
given explicitly in Table I. The same table also define the reduced matrix elements for 

TABLE I: Coupling constants. Here nfj = //p/6,/i]^ = 5/jp/6 with [ip = e/2m q . m q = 360 MeV 
is the quark mass which is determined here by including the one-loop corrections to fit the proton 
magnetic moment (see Table II). 



B' B S B > 


b(Tb'b) 


< B'\\S B 'b\ 


\B > IttB'B 
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JttNN 
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v/72/25/^tv 





^72/25/4 
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S A 


2\/l5 


fnNN 1 5 


/4 





isospin operators r for NN, T for AiV, T f for iVA, and T A for AA. 

Let us first calculate the mass correction terms T*° N and £ A that are given in Eqs. ()15|) - 
([TKjl. By using Eqs. (|HJ), (fbljl and (jSBI), we obtain in the rest frame of A" and A 

o r dk u f-wNN s 2 1 cr ■ kcr ■ kr ■ t \ F NtN (k) \ 2 

N ~ J (2tt)3 < m ^ m ^ l[[ m7T > 2E 7T (k) m N - E^k) - E N (k) 

+ { mj 2E n (k) m N - E n (k) - E A (k) J 1 m 
f dk ^ | / f n aa \ 2 1 5 A ■ fc5 A • fcT A ■ T A I F A , A (fc) | 2 , 

j (2vr)3 < msAmTA 11 m,, J 2£ 7r (fc) m A - ^(Jfe) - £7 A (fe) 1 ' 

(39) 

To perform the calculations, we need to define the form factor FB,B>{k) in Eq. ()33j) . 
To be consistent with the SL model, we here depart from the usual oscillator form and set 
F B:B ,(k) = (A 2 /(A 2 + fc 2 )) 2 with A = 650 (MeV/c) 2 for all nBB' vertices. Eqs. 
then lead to the following results 

Z% = T,° N (irN) + A) 

= -73.5MeV - 65AMeV = -139. MeV 
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£ A = £ a (ttA) 
= -76.6MeV 



(40) 



Here we indicate the intermediate state in each pion loop term. Note that irN intermediate 
state is excluded in the correction to the bare A mass, since its effect is already included in 
the rescattering term induced by the vertex interaction of Eq. (|29j) . The contribution 

of this rescattering to the mass shift is 

V res p f dk ,f*NAy 1 S^-kS-kT^ -T\F NA (k) | 2 

La " W (2^)3 < 1 ( ^7 j 2Ejk) m A - Ejk) - E N (k) 1 > 

= -46.2MeV (41) 

where P denotes taking the principal-value part of the integration. 

We now note that (Ho + + £ A ) of the effective Hamiltonian Eq. (|28|) defines the 
physical nucleon mass and the pole position W A = = 1232 MeV of the K-matrix of itN 
scattering in P33 channel. Thus, we have the following relations 

m N = m% + E° N , (42) 

and 

W A = m A + Y? A + £ r A es . (43) 

With the above results, the mass parameters m° N and m A associated with Hq of the 
effective Hamiltonian Eq. (|2*Hj) can then be determined 



m^ N = mN - J]% = W77MeV 

m ° A = W A — X A — Y>f s = 1355MeF 

The difference of the two bare masses are 

5° = m° A - m° N = 278MeV (44) 

The mass parameters m? N and m A obtained above can be considered as data for determining 
the parameters of a hadron structure model which 'exclude' the pion degree of freedom. 
Accordingly, one can assume that H of Eq. (|28|). which is defined by these two bare 
masses, can be identified with a model Hamiltonian defining the structure of the constituent 
quarks within the nucleon and A. 
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Most of the existing constituent quark model calculations [15j, |l6j, l^Ql, determine their 
parameters by fitting the mass difference Sm = m&(= 1232) — mjv(= 938.5) = 294 MeV, 
not by reproducing the absolute values of the masses of N and A. We note that this mass 
difference is not so different from that given in Eq. (J44j) . Thus we can identify H of Eq. 
(|28|) as the constituent quark model Hamiltonian with its eigenfunctions | N > and | A > 
consisting of three quarks. Accordingly, the current matrix element A defined by Eq. 
(|3*T|) can be identified with the prediction from the constituent quark models. 

We now turn to calculating the electromagnetic form factors. The constituent quark 
contribution is described by the current operator Jb,b' of Eq. (J3J). Within the considered 
SU(6) constituent quark model and in the second quantization notation of Eq. (jSJ), we can 
write 



J w 



B'B 



x (e£ + ^] P ^-8 BI , B + \n s B 6 B , tB + » V B , B TE, B }iS B ,, B x q)F B T B (q 



1 T z ~> i 

2 2 ' 2m B 

+ (h.c.)] (45) 

where F B m B ,(q 2 ) is an electromagnetic form factor, and the parameters fi B and ^i B are 
defined in Table I in terms of \xp = e/ (2m q ) with m q denoting the quark mass. In consistent 
with the SL model, the other two current operators in Eq. (jlj) are 

t7 J ( 27r ) 3 y/2E n (k) ^2E n (k') 

x {-iee^{k + k'))F^{q 2 ) + {Kc.)] (46) 



and 



J B > 



B',B,TT 



I d_P^dkdq 6{p _ q _ p ,_ k y, X ^^ bUp>Uk)bB{p 



W 12 * " " ' ^E~{k) 

_fnB',B rpj c \Tpem I „2\ 



x (e^^e t , 3 T B ,^ B S B , tB )F B 7 B7T (q 2 ) + (h.c)]. (47) 

Tfl 7,- 

The corresponding charge density operators are 

r dkdk'dg 1 1 iq . x 

^ h J ( 2 -) 3 v^M^v^O n {q) 

[S(k -q- k')al{k')a^{k){-iet^{E^k) + E^k'))) 
+ 5(k ~q + k')a^{k')a^{k){-iee ijZ {E^k) - E w (k'))) + (h.c.)], (48) 

PB>,B = Vfl/^^P(^ 
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TABLE II: Magnetic moment of nucleon in unit fj,^ 
'tree' Total with loop corrections 



proton 2.61 2.75 
neutron -1.74 -1.95 

In the above equations, the form factors F^p^q 2 ), F% m (q 2 ) and F B T Bn (q 2 ) should in prin- 
ciple be calculated from the associated hadron structure. This is a nontrivial task, as well 
recognized. For this exploratory study, we simply set all of these form factors equal to 
Fv(q 2 ) — V(l ~~ I 2 /My) 2 with M v = 0.76 GeV being the mass of vector meson. Obviously, 
this is the simplest prescription to maintain the gauge invariance. 

With the above definitions, we can evaluate loop corrections defined by Eqs. (J2U)-(|27J) 
by inserting appropriate intermediate states(illustrated in Fig. 7) and using Eq. (|T4"|) for Si 
and Eq. $T7$ for S 2 . 

To proceed, we need to first fix the quark mass m q which determine the current Jb'b- 
This is done by fitting the nucleon magnetic moments. The one-loop corrections (similar 
to what are shown in Fig. |7j) are included in the fit. We find that the nucleon magnetic 
moments can be reproduced very well if we set the quark mass as m q = 360 MeV. The 
results for the magnetic moments are shown in Table II. We see that the loop corrections 
are about 5% for proton and 10 % for neutron. It is important to note that the size of 
one-loop corrections depend heavily on the range A of the form factor Fb^b of Eq. (}33|) . 

We now turn to investigating the loop corrections on the — > A transition. Following 
the formulation presented in SL model, the 'jN — > A vertex function calculated in the A 
rest frame can be written in the following form 



< A I r 7A r^ A I q > 



E N (q)+m N 1 3(m A + m N ) 

(2tt)3/2^ 2E N (q) ^Am N (E N (q) + m N ) 3 

x [iG M {q 2 )S x q ■ e + G E {q 2 )(S ■ eer ■ q + S ■ qcr ■ e) 

, G c (q 2 ^ 



Sqaqeo], (50) 



where e = y An/137, q = (oj, q) is the photon four-momentum, and e = (e , e) is the photon 
polarization vector. The above definition allows us to calculate the multipole amplitudes of 
the 7Af — > A in terms of Gm (<? 2 ), Cs(g 2 ) and Gc(q 2 )- Explicitly, we havej^ 

A M (q 2 ) = [r 7 ^ A Wi = NG M (q 2 ) (51) 
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A E = [rV_> A ] fi2 = -NG E (q 2 ) (52) 
A c = [^ N -,a}c2 = N^G c (q 2 ), (53) 



with 



We first discuss 
determined in Refs. 



N _ e m A \q\ 1 

~ 2m N \ m N [l- g 2/( m7V + mA )2]i/2 \ > 

;he results at q 2 = photon point. The values of Am, Ae> and Ac 
are listed in Table III. These are the quantities we would like to 
interpret within the considered constituent quark model. Assuming that F 7 tv,a of Eq. (|3T)j) 
is what has been determined in the SL model, the values listed in Table III thus include 
the contribution not only from the quark-excitation term r° 7VA of Eq. (J31|) . which can be 
calculated from Eqs. and but also include the pion-loop contributions illustrated 
in Fig. [7| These loop contributions can be calculated by inserting appropriate intermediate 
states in the commutators of Eqs. (J2U)-(|21I)- As an example, we write down the expression 
for the mechanism Fig. 7b in the rest frame of A (p A = 0, p^ = —q) 



< A| - zSiJgyjiSiliV > = f dk J2 

J B=N,A 

< A|r°t |7r(fc), A >< A\J^ B \B >< 7r(fe),B|ro \N > 



(m A - E A (k) - E 7r (fc))( J E JV (qr) - E B (q + k) - E„(k)) 



Jdk £ 



B=N,A 

^a,a(0, -k, k)<A\A B \B> F wBtN (-k - q, -q, k) 



(55) 



(m A - E A (k) - E„(k))(E N (q) - E B (q + k) - E„(k)) 

An important point to note here is that the integrand in the loop-integration is independent 
of the collision energy W and has no singularity in the integration region < k < oo. Thus 
the included pion cloud effects are different from what were calculated from the SL model 
: [f 77r G' 7r Ar(W / )r 7r jv,A(W / )] which depends on the collision energy W in the ttN propagator 
G-kn(W). Qualitatively speaking, the one-loop contributions of Fig. are due to virtual 
pions which are part of the internal structure of N or A, while the SL model only accounts 
for the effects due to pions in scattering states which can reach the on-shell momentum 
asymptotically. 

The Q 2 = (Q 2 = —q 2 > 0) results from our complete calculations for all one-loop terms 
in Fig. [7|are presented in Table IV. In the first row, we list the values from r° Ar _ >A , which 
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is due to photon interactions with constituent quarks. As expected, the assumed spherical 
s-wave quark configurations do not have E2 and C2 transitions. In the same table we also 
list the contribution from each loop contribution illustrated in Fig. [7| The terms under 
'pion' and 'Seagull' are from Fig. 7a and 7c respectively. Fig. 7d gives the contribution 
'Normalization' which is the consequence of the appearance of the mass shifts Y>° N and £ A 
in the effective Hamiltonian Eq. ()3U|) and is naturally derived here by using the unitary 
transformation method. Fig. 7b contains the contributions due to spin transition and 
convection current. Their separate contributions are listed under 'Spin' and 'Convection' 
respectively. We note that Ac only has contribution from pion term because of the angular 
moment selection rule. 

We should emphasize here that the present calculations are based on a non-relativistic 
quark model and can only be compared qualitatively with the empirical values (table III) 
determined in the SL model. Thus we should not worry about their differences in absolute 
magnitudes. Rather, we focus on the relative importance between Am, Ae and Ac listed in 
Tables III and IV. 

The first interesting result in Table IV is that the total one-loop correction (Total - 
T^a) for the magnetic form factor Gm is only about 4 % of the 'bare' value r%- A , mainly 
due to the large cancellations between different contributions. In particular, the very large 
contribution from 'Normalization' of Fig. 7d plays a crucial role. We also see that the 
calculated Ae and Ac in Table IV are in opposite signs of the values listed in Table III. 
These results have the following implications. First, the T® N A values of Ae and Am in Table 
IV could be nonzero and negative such that the total values become the SL values listed in 
Table III. This can be the case if we assume that the quark wavefunctions of A and/or N 
could have a L = 2 d-state component. The other possibilities are that there could have 
multi-pion loop corrections and exchange current contribution of the quark electromagnetic 
current Q]. Within our formulation, some of these mechanisms can be derived from applying 
the third-order unitary transformation Us. A more detailed study along this line is clearly 
needed to make progress. 

The calculated Q 2 — dependence of the 'jN — > A form factors are displayed in Figs. |H1 
and 01 The dominant Ml transition is shown in Fig. |HJ The difference between the solid 
and dashed curves is due to the one-loop corrections. It is very weak and only visible at very 
low Q 2 . On the other hand, the pion cloud effects due to scattering states (dashed curve) 



17 



TABLE III: 'bare' helicity amplitudes of the SL model[l,2]. Unit is lO^GeU^ 1 / 2 . 



A M 173.3 
A E -2.3 
A c -2.2 

TABLE IV: 'bare' helicity amplitudes in quark model with loop correction. Unit is W^GeV^ 1 ^ 2 . 

SL is the result from Ref.[l,2] 

rJyiVA Pion Spin Convection Seagull Normalization Total 
A M 204.9 23.7 24.2 -3.3 -35.6 213.8 

A E 2.8 -0.1 0.1 0.9 3.6 

A c 1.1 0.0 1.1 

is very large, As discussed in detail in Ref. jjj, this finding explains why the conventional 
constituent quark model predictions disagree with the empirical value of the magnetic Ml 
transition of — ► A. The present result for one-loop corrections do not change that 
conclusion. 

The situation for Ae(Q 2 ) and Ac(Q 2 ) is quit different. Here we do not have contribution 
from quark excitation term r° JVA because of the assumed L=0 wavef unctions. We see 
that the calculated one-loop corrections (solid curves) are comparable in magnitudes to the 
pion cloud effects due to scattering state (dashed curves) calculated in SL model. More 
importantly, they have very different Q 2 -dependence and are opposite in signs. As seen in 
Eq. (|30p. the solid curves must be interpreted as part of the bare form factors determined 
phenomenologically in SL model. Namely, the form factors obtained from subtracting the 
solid curves from SL model's bare form factors are the contribution from quark excitation. 
This will be an important information for testing various hadron structure calculations. 
However, such information can not be realistically extracted here because of the simplicity 
of the model employed. 
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FIG. 8: Am{Q 2 )- Dot-dashed is from quark excitation A , the solid curve is the sum of F® N A 
and fI~^°£ p . The dashed curve is from pion scattering calculated in SL model. 
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FIG. 9: Left:AE{Q 2 ), Right:Ac{Q 2 )- Solid curves are from one-loop corrections, while the dashed 
curves are from pion scattering calculated in SL model. 

IV. ONE-LOOP CORRECTIONS ON NON-RESONANT 7 A irN 



One of the difficulties the SL model has in describing the data is from the non-resonant 
amplitude, fn this section, we would like to explore whether this can be improved by 
including the one-loop corrections vl~ loop of Eq. (}32|) . As a start, we will focus on the near 
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threshold region and consider only the the E Q+ amplitude. A complete calculation of v l~ loop 
for all partial waves up to A resonance energy is much more involved and will be explored 
elsewhere. 

First, we point out that the SL model failed to describe the near threshold 7p — > 7r°p 
data. For example, at = 145 MeV the SL model gives (after taking into account the 
effects due to the mass difference between 7r° and ^ ) 

E 0+ (U5MeV) = -2A7(Born) + 2.31(Rescattering) = -0.15[l(T 3 /"v+] (56) 

where Born is from the non-resonant production operator v^ e constructed in SL model, 
Rescattering include the effects due to final ttN interaction. The empirical value is E^ (145 
MeV) ~ -1.50. In getting the above result, we find that the main contributions to the Born 
term are from the nucleon-direct and nucleon-exchange diagrams, while the rescattering 
term is mainly from pion-pole and contact interaction through , -f+p^7i + + n—>-7r +p 
charge-exchange process. We also find that the s-wave charge exchange pion rescattering 
is dominated by the p-exchange n — N potential and the Born approximation t^N ~ v-kn 
is accurate. Furthermore the short range approximation of p-exchange potential (l/{m 2 + 
{Pn — Pn) 2 ) ~ l/ m p) 1S accurate within 10% in determining the rescattering effects in the 
considered near threshold energy region. With these considerations, the one-loop corrections 
near threshold can be calculated with the following much simplified Hamiltonian 

Hi = JL ^-^Nlhlnd IJl (j) 1 , ■ fip N + XipN^TipN ■ 07T x . (57) 

777-7,- 

Here the second term is a contact interaction with the strength determined from the p- 
exchange coupling constants : A = g p7T7T g p NN/ (2m 2 p ). By minimum substitution, the second 
term of Eq. (|57jl will generate a interaction current 

3N,Nwn = e\[i> N ^Ti) N X $ v ] X ^ . (58) 

It induces an electromagnetic contact interaction involving two pions. To maintain the gauge 
invariance within the model defined by the simplified interaction Hamiltonian Eq. ([57)1 . this 
current is included in the calculation along with the currents J n , Jb',b and Jb',Btt given in 
Eqs. (|46|) - (j48|) and illustrated in Fig. |21 All coupling constants and vertex form factors are 
taken from SL model. 

With the above simplified model, we first re-calculate the rescattering contributions, 
~ v t J^ e GnN(W)v n N to the E + amplitude for 7p — > The results at E 1 = 145 MeV are 
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listed in Table V. It is instructive to note here that the calculated rescattering contribution 
involves cancellation between the terms (d) and (e). The total rescattering value 2.36 is very 
close to the value 2.31 of the rescattering term in Eq. of the SL model. This justifies 
the use of the simplified model defined by Eqs. (J57|) and (J58j) . 

TABLE V: Rescattering contributions to the E + amplitude of ~yp — ► ir°p at 145 MeV, calculated 
from mechanisms (b)-(e) illustrated in Fig. using the model defined by Eqs. (|57 |) -(|58 |) . 

Diagram (b) (c) (d) (e) sum 
-0.074 -0.685 -1.966 5.087 2.36 



/I 

(a) 



(b) 



(c) 



(d) 



FIG. 10: Subset of loop corrections on the — > ttN transition amplitude. 



(a) 



/I 



(b) 



(c) 



(d) 



(e) 



A 



A 



(f) 



(g) 



(h) 



(j) 



FIG. 11: Subset of Loop corrections on the jN —* 7T./V transition amplitude. 



The one-loop corrections can be calculated from Eq. (j25j) by inserting appropriate in- 
termediate states. The resulting amplitudes are illustrated in Figs. El and EI Note that 
diagrams in Figs. El and E] are not time-ordered diagrams. Rather they just illustrate the 
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structure of the matrix element of each term in t>^~' oop . As seen from Eq. (jl4j) and Eq. ()17j). 
the loop integrations for all processes in Figs. El and EH will involve energy — independent 
propagators associated with these two operators. Thus, although the diagrams in Fig. [TU] 
look similar to the rescattering terms, but they are energy — independent matrix elements. 
The calculations for these loops are tedious but straightforward, and will not be elaborated 
here. 

TABLE VI: one-loop contributions to the E + amplitude of "fp — > 7r°p at 145 MeV, calculated 
from mechanisms illustrated in Fig. 11 using the model defined by Eqs.(54). 







Diagram 




Diagram 




Fig. 


. 10 


(a) 


-0.079 


(c) 


-0.024 






(b) 


-0.090 


(d) 


0.475 


Fig. 


. 11 


(a) 


0.157 


(f) 


-0.418 






(b) 


-1.192 


(g) 


0.00 






(c) 


0.875 


(h) 


0.00 






(d) 


-0.085 


(i) 


-0.696 






(e) 


-1.011 


0) 


0.699 



Sum = -1.39 

Our results at = 145 MeV for each of the one-loop corrections shown in Figs. 10 and 
11 are listed in Table VI. The results listed in Tables V and VI lead to 

E 0+ (U5MeV) = -2A7{Born) + 2.32(Rescattering) - 1.39{Loop) = -1.54 (59) 

This reproduces the empirical value Eq^ (145 MeV) ~ -1.50. The calculated effect of the 
one-loop corrections for E 0+ in the near threshold energy region is shown in in Fig. IT21 
Clearly, the one-loop corrections drastically reduce the magnitudes and bring the results to 
agree with the empirical values. The kinks due to the cups effect are reproduce well in our 
calculations. 

In Fig. we show that the one-loop corrections on the Eq+ amplitude can change 
significantly the calculated angular distributions to better agree with the data. To see the 
full one-loop correction effects, we need to also calculate other multipole amplitudes. This 
along with the results for the A region will be explored elsewhere. 
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FIG. 13: 7p — > tt°p at 146, 150, 160, and 170 MeV. Solid curves are obtained when the one-loop 
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V. SUMMARY AND OUTLOOK 



In this paper, we have applied the unitary transformation to derive the leading-order 
corrections on the effective Hamiltonian of the SL model for electromagnetic pion production 
reactions. We have investigated the one- loop corrections on the masses of N and A, the 
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7^V — > A vertex, and the non-resonant pion production operators. Qualitatively speaking, 
the derived one-loop corrections are due to the virtual pions which are part of the internal 
structure of iV or A, while the pion cloud effects generated within the SL model or the 
other dynamical models, such as the Dubna-Mainz- Taipei (DMT) model j^, only account 
for the effects due to pions in the scattering states which can reach the on-shell momentum 
asymptotically. 

With the one-loop corrections included in determining the mass parameters, we find that 
the free Hamiltonian of the model can be identified with the conventional constituent quark 
model. We then proceed to apply such a constituent quark model to calculate the one-loop 
corrections on the 7iV — > A transition form factors. It is found that the one-loop corrections 
on the magnetic Ml transition is very small. Our results further establish the conclusion 
reached by the SL model that the large discrepancy between the conventional constituent 
quark model predictions and the empirical values are due to the pion cloud effects associated 
with the pions in scattering states. 

The calculated one-loop contributions to the electric E2 (Ae) and Coulomb C2 (Ac) 
form factors of the 'yN — > A transition are found to be in opposite signs of that due to pion 
cloud associated with the scattering states. One possible implications of this result is that 
the extracted empirical values of SL model could be largely due to the nonspherical L = 2 
intrinsic quark excitations which could lead to nonzero and negative contributions to Ae and 
Ac- On the other hand, there could have higher-order exchange current contributions which 
are not included in this work, but must be also calculated for a complete understanding of 
the empirical values of SL model. Clearly more works are highly desirable. 

We have also found that the one-loop corrections on the non-resonant pion production 
operator can resolve the difficulty the SL model encountered in reproducing the empirical 
E 0+ amplitude of near threshold 7r° photoproduction. It will be worthwhile to further 
extend this work to calculate these one-loop corrections for higher partial waves. Some of 
the discrepancies between the SL model and the data in the A excitation could be removed 
by including these corrections. Our effort in this direction will be reported elsewhere. 

To end, we emphasize that the SL model is obtained from keeping only the lowest order 
terms of a formulation within which the higher order terms can be rigorously derived. At- 
tempts to fit the data by adjusting the current SL model are not justified theoretically. The 
most important task to improve the SL model is to include these corrections order by order 
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until the convergence of the predictions has achieved. In this work we have taken a very first 
step in this direction. Undoubtly, much more works are needed to complete a consistent 
dynamical model of electromagnetic pion production reactions. 
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